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We study the one-dimensional S = 1/2 Heisenberg model with a uniform and a staggered 
magnetic fields, using the dynamical density-matrix renormalization group (DDMRG) 
technique. The DDMRG enables us to investigate the dynamical properties of chain with 
lengths up to a few hundreds, and the results are numerically exact in the same sense as 
'exact diagonalization' results are. Thus, we can analyze the low-energy spectrum almost 
in the thermodynamic limit. In this work, we calculate the dynamical spin structure 
factor and demonstrate the performance of the DDMRG method applying the open-end 
boundary condition as well as the periodic boundary condition. 
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1. Introduction 

In one-dimensional electron systems, the conformal field theory has succeeded in 
the description of the low-energy physics of the Tomonaga-Luttinger liquid. Recent 
experimental progress, such as angle resolved photoemission and neutron scattering, 
demands the understanding of the elementary excitation pictures of the dynamics 
not only in the low-energy and long- wave length limit, but over the wide range of 
frequency and momentum. 

In no magnetization case, Miiller et al^ proposed the approximate conjecture 
of the spin structure factor S{q,uj) in the one-dimensional (ID) S ~ 1/2 antifer- 
romagnetic Heisenberg model with nearest neighbor interaction. By usage of the 
mathematical works, the exact expression for two-spinon contributions to S{q, to) 
in the thermodynamic limit was found. ^ The exact result and Miiller ansatz for 
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S{q,w) have the same singularity at the lower spectral boundary, which is called 
the des Cloizeaux- Pearson mode. Remarkable progress has been made on the Bethe 
Ansatz (BA),^ which enables us to numerically calculate the S{q, uj) for the Heisen- 
berg model in a magnetic field for relatively large system size.^ 

Recently, a dynamical density-matrix renormalization group (DDMRG) method 
has been developed for calculating dynamical correlation functions at zero tem- 
perature in quantum lattice models. It is based on a variational formulation of 
the correction vector technique. ^'^ This method is an extension of the standard 
DMRG method^ which is a powerful numerical technique for a variety of ID sys- 
tems, whereby we can obtain very accurate ground state and low-lying excited 
states. So far, it has been shown that the momentum-dependent dynamical quan- 
tities can be calculated with good resolution in finite open Hubbard chains by the 
DDMRG method.^ We could also expect comparable performance for the same 
model in magnetic fields. In this work, we apply the DDMRG method to the ID 
S = 1/2 Heisenberg model with an uniform and a staggered magnetic fields and 
calculate the spin structure factor S{q,u}). The accuracy of the DDMRG method 
is checked by comparing to some analytic solutions. In addition, the difference of 
performance by the boundary conditions is demonstrated. We note that the Heisen- 
berg model can be dealt with the DMRG technique much easier than the Hubbard 
model. 

This paper is organized as follows. In the next section, we introduce the model 
and the spin structure factor. Section [3] is devoted to show the brief outline of 
the DDMRG method. In Section [H we demonstrate the DDMRG calculations of 
the spin structure factor for various magnetic fields. The performance with open 
boundary condition (OBC) as well as with periodic boundary condition (PBC) is 
examined in comparison with the exact solutions. Some remarks are given in the 
final section. 

2. Model 

We consider the ID S = 1/2 Heisenberg model in magnetic fields: 

H^jY^S^- S,+i -HY^St + hY{-^ySt, (1) 

i i i 

where J(> 0) is an antiferromagnetic exchange interaction between neighboring 
sites. Henceforth we take J as the unit of energy. H and h are a uniform and a 
staggered magnetic fields along the z-direction, respectively. 
The dynamical spin structure factor is defined as 

n>0 

where {— 5*"^^) is the Fourier transform of the spin operator S'f at site i, and En 
and IV'n) are, respectively, the n-th eigenenergy and eigenstate of the system (with 
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the ground state denoted by n = 0). Note that the definition of the momentum- 
dependent spin operator depends on choosing the boundary conditions. The 
definitions are given in Sectional 



3. DDMRG method 

In this section, the DDMRG method is briefly summarized for calculating the 
dynamical spin structure factor Thus, we are interested in the momentum- 
dependent spectral functions, 

Aiq,Lu) = lunJ-lm{MS%-^ ^— -^,^|^o). (3) 

The dynamical spin structure factor is then obtained as 

S{q,i^) = --lTsiA{q,u) 

TT 

= limIm(.^o|5^,-^ "^^.^ , (4) 

To compute the dynamic quantity such as Eq. ([3]), we use the DDMRG method. 
This approach is based on a variational principle. One can easily show that for 77 > 
and fixed frequency lo the minimum of the functional 

with respect to all quantum states IV') is 

2 

W'(V'mi„) = {MS%7^ ^\r. — (6) 

The functional minimum is related to the convolution of the dynamical spin struc- 
ture factor with a Lorentz distribution of width 77 by 

M^(^min)--^?7^(g,C^). (7) 

The DDMRG method consists essentially of minimizing the functional ([5]) nu- 
merically using the standard DMRG algorithm. Thus the DDMRG provides the 
dynamical spin structure factor S^{q,oj) for a finite broadening 77. The full spec- 
trum ([4]) convolved with the Lorentz distribution is given as 

/oo 
Auj'S{q,uo')— % (8) 

The real part of Eq. Q can be calculated with no additional computational cost 
but is generally less accurate. The necessary broadening of spectral functions in 
DDMRG calculations is actually very useful for studying continuous spectra or for 
doing a finite-size scaling analysis.'' 

If one would like to obtain the spectrum in the limits of L — > c» and 77 = 0, 
it can be done by carrying out a deconvolution of the DDMRG data.^° In theory, 
the deconvolution scheme corresponds to solving ([8]) for obtaining S(^q,uj), where 
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Fig. 1. 2-component of the total-spin quantum number as a function of H (a), and as a function 
of h for various fixed H (b). The solid line denotes the exact behavior for H — > 0. 

a set of S^QjCli) on the left-hand side is calculated with the DDMRG method. 
We also know that the broadened spectrum of the system on an infinite lattice 
(L — > oo) is usually almost identical to the spectrum on a finite lattice (N < oo) 
if ?7 > c/L (the constant c is comparable to the width of the spectrum). Therefore, 
assuming that the DDMRG data S^'{q,u!) describes the broadened spectrum for 
N oo, one can solve (|S]) approximately under the condition that S{q,oj) is the 
convolved 'exact' spin structure factor. We note that S^{q,uj) must be a continuous 
and relatively smooth function in order to obtain quantitatively accurate spectrum 
after deconvolution. To this end, it is required to choose rj smaller than the width 
of the spectrum in 77 = 0.^ 

4. Results 

4.1. Field-induced magnetization 

The z-componcnt of the total-spin 5**°* is now a good quantum number for the 
DMRG calculation. Therefore, it is very efficient to classify the Hamiltonian by 
5**°'. Let us first study 5**°* of the ground state in the system with magnetic fields 
before starting the dynamical calculations. In Fig.[lja), we show S'*°* of the ground 
state as a function of the uniform magnetic field H. The calculations are carried 
out with the DMRG method in a chain with i = 40 sites under the PBC. Note that 
5**°* is equivalent to the magnetization of the system via the relation M = 31°*" /L. 

For small H, the magnetization is induced proportionately with the uniform 
field, M = {2/tt'^)H}^ The derivative ^ becomes larger with increasing H since 
the antiferromagnetic spin fluctuations are rapidly weakened with approaching to 
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Fig. 2. Spin structure factor S+~{q,uj) (a) and ^^^(g.aj) (b) for L = AO, H = 1.537, and t) = 0.1 
applying the PBC. Dotted curves show the BA solution with r} = O.l for comparison: up to two- 
spinon and two-antispinon contributions (a) and up to two-spinon and one-antispinon contributions 
(b). 



the fully-polarized phase. The system is fully polarized at a critical uniform field 
He = 2 where the spin susceptibility diverges, i.e., oo. We also show the 

dependence of 5*'°* on the staggered magnetic field h for several kinds of H values in 
Fig.[ljb). We can see that the magnetization is rapidly suppressed by the staggered 
field and it reaches to zero at a critical value he- The critical staggered field is 
roughly scaled as he ^ H/2. 



4.2. Dynamical spin structure factor 
4.2.1. Periodic boundary condition 

Normally, the OBC is applied to the (D)DMRG calculations in order to achieve an 
accurate treatment; namely, we can study systems with the OBC on much larger 
lattices than systems with the PBC, with keeping a sufficient accuracy. Especially 
in interacting fermion systems such as the Hubbard and the t—J models, the OBC 
seems to be rather essential to adequately truncate their large Hilbert space. How- 
ever, a balk system is generally better described using the PBC, and moreover 
the exact definition of the momentum-dependent operators with the OBC is quite 
difficult. 

Here, an idea will naturally arise as follows: it may be that the PBC can be 
applicable in a spin system where the Hilbert space is much more reduced than that 
in the fermion systems. Motivated by this idea, we check the DDMRG performance 
with the PBC in the ID 5 = 1/2 Heisenberg model. We now start to calculate 
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Fig. 3. Average discarded woigiit as a function of system size L for the truncation number 
m = 100, 200, and 400 in the calculation of S'^ (0, uj). Dotted lines are guides for eyes. 

the dynamical spin structure factor With the PBC, the spin operators are 
defined by 

with momentum q = 2ttz/L for integers — i/2 < z < L/2. 

Using BA sohitions on a finite lattice, the spectral weight are available in the case 
with the uniform magnetic field. "^'^ Comparison with the BA solution will provide an 
opportunity to test the performance of the DDMRG method. In Fig. [2l we compare 
the DDMRG spectra with the line shapes based on the BA solutions in iJ = 1.537, 
where the system size is L = 40 sites and the broadening is taken as 77 = 0.1 for 
both of the methods. We can see excellent agreements in the low-energy excitations, 
whereas the BA spectra deviate from the DDMRG data in the high-frequency range. 
It means that the DDMRG method can take into account higher-order terms than 
the BA treatment. If the discarded weight is sufficiently small in the DMRG 
truncation procedure, e.g., Wd < 10~^, the DDMRG results for finite systems are 
numerically exact in the same sense as 'exact diagonalization' results are. In the 
present calculations, we keep up to m = 300 density-matrix eigenstates and the 
typical discarded weight is less than 5 x 10~^. Thus, we are confident that the 
DDMRG method using the PBC works well in the iS = 1/2 Heisenberg chain with 
at least up to several tens of sites. 

Let us now mention how large system can be dealt with using the DDMRG 
method. The performance depends mainly on the number of the density-matrix 
eigenstates truncated m. As the truncation number m increases, the accuracy of 
calculation is improved but more CPU time and disk space are required. Practi- 
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cally, we could keep a several hundreds density-matrix eigenstates with the present 
typical cluster machine, e.g., Pentium 4 3.2GHz or Opteron 252 2.6GHz. In Fig. [31 
we show the logarithm of average discarded weight Wd as a function of system size 
L for m = 100, 200, and 400 in the calculation of (0,cj). Naturally, Wd de- 
creases as m increases at a fixed system size. On the other hand, for a fixed m, Wd 
increases rapidly with increasing L for smaller systems and is almost constant for 
larger systems. Empirically, we can no longer expect a good DMRG performance 
in this 'constant' region. For example, we estimate L = 40, 60, and 80 as the opti- 
mum system sizes for m = 100, 200, and 400, respectively. Therefore, the DDMRG 
calculation with about a hundred lattice sites would be adequately possible in the 
standard S* = 1/2 Heisenberg chain using the PBC if we take m « 600. 

We then add the staggered magnetic field to the system under a fixed uniform 
magnetic field. In Fig. IH we show the spin structure factors S°'^{q,u}) for af} = 
+-, — h, zz at = 0, 0.08, 0.16, 0.24, and 0.31 with fixing H = 1.6, calculated with 
the DDMRG method for i = 40 and rj = 0.05 applying the PBC. The high-energy 
structures, which cannot be obtained by the conformal field theory, are clearly seen. 

For h — [see also Fig. [2], the feature of S{q,u!) is very similar to that of 
Haldane-Shastry model (l/r^ interaction). ^^"^'^ In the lower energy edge of the small 
q {0 < q < 2ttSI°^/L) in the (9,^), the one antispinon (magnon) contributes as 
a (5-function peak. At g = 0, the exact form is known, (0, uj) = (47rS'*°*/i)^(w — 
H) , which presents a resonant mode of the magnetization induced by the uniform 
magnetic field0 The spectra are changed gradually with increasing h and we can 
see some predominant features for h — 0.32 as the following: 

• The low-energy structures of ^\q,uj) around q = arise. They are associ- 
ated with increase of the weight of unfluctuating antiferromagnetic configuration 
in the ground state. 

• The one-magnon peaks of (q, lo) in the small momentum region are smeared 
and they will disappear when h > he- 

• In the large momentum {q ~ tt) region of {q,'-^), the continuum structure 
when h ~ turn to a few discrete peaks due to the decrease of spin fluctuations. 

• The band structure of S'^^{q,uj) seems to be split into two small-dispersive bands 
around w < 0.5 and oj ^ H. 

4.2.2. Open-end boundary condition 

We now turn to the calculation with the OBC. A DDMRG resuh with the OBC 
is less exact than that with the PBC for the same system size. However, the OBC 
enable us to carry out a calculation for much larger systems, and then we may obtain 



''The one-antispinon and the two-spinon plus two-antispinon contribution to (g, ui) in Haldane- 
Shastry model can be expressed in terms of the analytic expression of the advanced Green function 
for the spinless Sutherland model with coupling parameter A = 2.^^^ Other components, S ^(q,u}) 
and the two-spinon plus one-antispinon contribution to S^^{q,Lj) for Haldane-Shastry model also 
can be expressed by the correlation function of the spinless Sutherland model. ^'^'^^ 
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Fig. 4. Spin structure factor 5^ (9,1^) (upper), S ^{q,uj) (middle), and 5^^(9,0)) (bottom) 
calculated with the DDMRG method for L = 40 and 77 = 0.05 applying the PBC. The staggered 
magnetic field h is varied with fixing the uniform magnetic field at H = 1.6. 



the dynamical quantities in the limits of L ^ 00 and 77 ^ through additional 
techniques such as the finite-size scaling or the deconvolution techniques. Hence, it 
would be also useful to check a performance with the OBC. 

When the OBC is applied, we usually use the eigenstates of the particle-in-a-box 
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Fig. 5. Spin structure factor S^^{q,ui) in a chain with L = 192 sites for q = (192/193)7r ~ 7r (a) 
and q = {97/193)7r ~ it/2 (b), calculated with the DDMRG method for rj = 0.05 (dotted lines) 
applying the OBC and then deconvoluted (bold lines). The thin lines show the analytic solutions 
in the thermodynamic limit. ^ 



problem to define the operators 

s'^ = J^-Y,sin{qi)sr, (10) 

with quasi-momcntum g = 7rz/(L + 1) for integers 1 < z < L. This definition of 5^ 
should be equivalent to that for the PBC in the thermodynamic limit, but the agree- 
ment may not necessarily good in finite systems. In Fig. [5l we show the DDMRG 
results of S^^{q,u!) without magnetic field for L = 192 and rj = 0.05 with OBC, as 
well as the exact line shape of the two-spinon excitation contribution for L — )■ oo."^ In 
order to compare the finite-size DDMRG spectrum with the exact solution without 
broadening, we need to eliminate the broadening 77 from the 'convolved' DDMRG 
spectrum. 

An approach for obtaining the spectrum in the 77 = limit is the deconvolution 
of the DDMRG spectrum. Theoretically, a deconvolution means solving ^ for 
S{q, ijj) using the DDMRG data in the left-hand side. If it was possible to do this 
calculation exactly, one would obtain the discrete spectrum on a finite lattice of L 
sites. Nevertheless, the broadened spectrum on a infinite lattice {L — > 00) is usually 
almost identical to the discretized DDMRG spectrum (L < 00) under the condition 
77 > c/L. Since the width of S^^{q,Lo) is always less than c — 4 according to the 
exact solution, our choice 77 = 0.05 indeed satisfies the condition. 

At (7 = (192/193)7r, the deconvoluted DDMRG data seems to agree well with 
the exact solution in the presented scale [see Fig. [5] (a)]. Actually, the deconvo- 
luted DDMRG data is a shade thinner than the exact solution. It might come from 
the difference of momentum taken in the calculations, i.e., q = (192/193)7r in the 
DDMRG and tt in the exact solutions. Nevertheless, it can be rather hard to find 
the existence of l/cj-divergence for the lower edge uj = expected from the exact 
solution. If the width of a peak for 77 = is smaller than 77, the we cannot obtain it 
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Fig. 6. (a) Spin structure factor S^^{q,uj) at q = {L — 1/L)n ~ tt for various system lengths L 
and r] = 9.6/L. (b) Extrapolation of Smax to — > 0. (c) Extrapolation of -Emax to r; ^ 0. 



accurately by the deconvolution technique. Hence, we could not see the very sharp 
structure of S^^{tt,uj) at w < 0.05 correctly. On the other hand, the agreement is 
not very good a.t q — (97/193)7r ~ tt [see Fig. [5] (b)]. Occasionally, quantitative 
estimation would be somewhat problematic. We find that the accuracy of spectrum 
is usually the worst around q = 7r/2 for using the OBC. However, we suggest that 
the definition of operators with the quasi-momentum gives a good approximation 
at least qualitatively. 

The other approach is the finite-size scaling analysis of the DDMRG data.^'^ In 
Fig.[6](a), we show the spin structure factor S^^{q,uj) a,t q = (192/193)7r for several 
system sizes L. We keep a relation r/L = 9.6 for systematic extrapolation. We can 
clearly see the convergence of the finite-size spectra toward the exact spectrum 
as 77 decreases. We also can study the lower-edge behavior more quantitatively by 
a scaling analysis of the maximal value in the DDMRG spectrum S^^{q,uj). For 
instance, the height of the low-energy maximum logj^g Smax is scaled linearly as a 
function of log^g r] [see Fig. |6] (b)] and the slope is —1. It means that the spectrum 
diverges as i]"^ for rj —> 0. Moreover, the position of the low-energy maximum, 
EmuK, approaches to as 77 goes to 0. We can then confirm that S^^{tt,uj) has a 
singularity with exponent —1 at w = in the thermodynamic limit L 00. This 
is consistent with the exact result which shows the l/w-divergence of S^^{tt,u!) at 
the lower edge u; — 0. 



5. Remarks and Outlook 

In this work, we demonstrate successful application of both the OBC and the PBC 
to the DDMRG calculation in the ID S = 1/2 Heisenberg model. Each of the 
boundary conditions has advantages and disadvantages, and then we need to make 
proper use of the boundary conditions depending on the situation. 
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(i) With the PBC, the system size must be restricted practically up to about a 
hundred but the result is numerically exact. Therefore, the PBC is suited to 
study 

- quantitative estimation of the whole spectrum, and 

- accurate picture of complex dispersive structure. 

(ii) With the OBC, we can set the system size to be two hundreds or more, and 
thus we can easily obtain the spectrum in the thermodynamic limit via the finite- 
size scaling or the deconvolution techniques. However, quantitative accuracy is 
occasionally missing on a finite lattice. Therefore, the OBC is suited to study 

- the behavior and position of a singularity in the spectrum, and 

- rough picture of simple dispersive structure. 

Possibly, the finite-size scaling or the deconvolution techniques will be applicable to 
the DDMRG data with the PBC if the spectrum had relatively simple form. 

Finally, we add that the method can easily be extended to other systems: for 
example, 

- 5=1/2 Heisenbcrg ladder systems, 

- more realistic models which include the transverse staggered magnetic field or the 
Dzyaloshinsky-Moriya interaction, and 

- spin systems where the total spin is more than S = 1/2, etc. 
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